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We have considered the conductivity properties of a two dimensional elec- 
tron gas (2DEG) in two different kinds of inhomogeneous magnetic fields, i.e. 



> 

f^ , a disordered distribution of magnetic flux vortices, and a periodic array of 



magnetic flux vortices. The work falls in two parts. In the first part we show 
how the phase shifts for an electron scattering on an isolated vortex, can 
be calculated analytically, and related to the transport properties through 



'^ ' a force balance equation. In the second part we present numerical results 



for the Hall conductivity of the 2DEG in a periodic array of flux vortices. 
We find characteristic peaks in the Hall conductance, when plotted against 



H ' the filling fraction. It is argued that the peaks can be interpreted in terms 

of "topological charge" piling up across local and global gaps in the energy 

spectrum. 
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I. INTRODUCTION 

Over the last decade the two dimensional electron gas (2DEG) have been exposed to a 
wide range of physical experiments, in which the electrons have been perturbed by differ- 
ent configurations of electrostatic potentials, with or without a homogeneous perpendicular 
magnetic field. These experiments have shown new kinds of oscillations in the magnetocon- 
ductivity, with a periodicity not given by the geometry of the Fermi surface, as is the case 
with the Shubnikov-de Haas oscillations, but given by the interaction of the two length scales 
given respectively by the magnetic length, and by the spatial structure of the potential, e.g. 
the Weiss oscillationsEJ. More recently, there have been increasing interest in systems where 
the 2DEG is exposed to an inhomogeneous perpendicular magnetic field. In such systems 
the inhomogeneities in the magnetic field acts as perturbations of the 2DEG, relative to the 
homogeneous magnetic field, where the band structure consists of the dispersionless Landau 
bands. The inhomogeneous magnetic field appears in the Hamiltonian in the form of a non 
trivial vector potential. In the case of a periodic variation in the magnetic field, it is possible 
to construct a periodic vector potential, if and only if the flux through the unit cell of the 
field is equal to a rational number, when measured in units of the flux quantum 0o = h/e. 
Under these special circumstances Bloch states can be used as a basis for the calculation of 
response properties of the electron gas. 

In this paper we have considered a special class of spatially varying magnetic fields 
which consists of flux vortices, that are either distributed at random or placed in a regular 
lattice structure. A system consisting of a 2DEG penetrated by a random distribution of 
magnetic flux vortices, have been experimentally realized by Geini et al.Bd. They made 
a sandwich construction of a GaAs/GaAlAs sample with a 2DEG at the interface, and 
a type II superconducting lead film (electrically disconnected from the 2DEG). When the 
system was placed in an external magnetic field, and cooled below the transition temperature 
of the film, the magnetic field penetrated the film, and thereby the 2DEG, in the form of 
Abrikosov vortices. When the external magnetic field is weak, below lOOG, the vortices will 



be well separated, and the 2DEG therefore sees a very inhomogeneous magnetic field. In the 
experiments conducted by Geim et al. the fiux pinning in the film was strong, resulting in a 
disordered distribution of flux vortices. This is the physical situation which we investigate in 
Sec. |1| below. In very clean films of type II superconducting material, the fiux vortices will 
order in a periodic array, i.e. an Abrikosov lattice, and thereby create a periodic magnetic 
field at the 2DEG. This is the situation which we analyse in Sec. |IITi 

Several authors have investigated the transport properties of 2DEG's in different kinds of 
inhomogeneous magnetic fields. Peeters and Vasilopoulosa have made a theoretical study of 
the magnetoconductivity in a 2DEG in the presence of a magnetic field, which was modulated 
weakly and periodically along one direction. They found large oscillations in the longitudinal 
resistivity as a function of the applied magnetic field strength. These oscillations are due to 
the interference between the two length scales given respectively by the period of the lateral 
variation of the magnetic field, and by the magnetic length corresponding to the average 
background field. The oscillations are reminiscent of the Weiss oscillations, but have a higher 
amplitude and a shifted phase, relative to the magnetoresistance oscillations induced by the 
periodic electrostatic potential. 

The problem of how the transport properties of the 2DEG is modified by the presence 
of a random distribution of fiux vortices, have been treated earlier by A. V. KhaetskiiH, and 
also by Brey and FertigH. The approach used by these authors are closely related to the 
one we have presented in Sec. 0, i.e. based on the scattering theory, and the Boltzmann 
transport equation. Khaetskii considered the scattering in certain limiting cases, including 
the semiclassical, and Brey and Fertig calculated the transport properties numerically. 



In Sec. |TT| we will address the "paradox" of how the Hall effect can disappear in the 
following situation: We imagine a 2DEG in a regular 2D-lattice of fiux vortices, with the 
magnetic field from a single vortex exponentially damped with an exponential length ^, 
in units of the lattice spacing. We take the total fiux from a single vortex to be 0o/2, 
as is the case when the vortices come from a superconductor. When ^ ^ 1, the field is 
homogeneous and the Hall conductivity is an = pe^ /h, where p is the number of filled 



bands. In the other hmit i.e. when ^ — i> 0, the time reversal symmetry is restored, and the 
Hall conductivity vanish. The fact that the system has time reversal symmetry when the 
vortices are infinitely thin, can be seen by subtracting a Dirac string carrying one quantum 
of magnetic fiux (pQ, from each fiux vortex. The introduction of the Dirac strings can not 
change any physical quantities, and the procedure therefore establishes that the system, 
with infinitely thin vortices with fiux +0o/2, is equivalent with the system with reversed 
fiux — 0o/2, through each vortex. The paradoxical situation arises because it is known from 
general argumentg3, that the contribution to the total Hall conductivity from a single filled 
nondegenerate band is a topological invariant, and therefore cannot change gradually. The 
situation is even more clear cut if we imagine a periodic array of exponential fiux vortices 
carrying one fiux quantum 0o each. Then the limit ^ = 00 corresponds to a homogeneous 
magnetic field, while the opposite limit ^ = coresponds to free particles. Incidentally this 
scheme can be used to establish an interpolating path between the multi-fractal structure 
known as Hofstadters butterfiyH which is a plot of the allowed energy levels for electrons on 
a lattice in a homogeneous magnetic field, and the corresponding plot for lattice-electrons 
in no field, which is completely smooth. 

The plan of the paper is as follows. In Sec. y we calculate the longitudinal and transverse 
resistivities of the 2DEG, in the limit where the mean distance between the vortices is so 
long that correlation effects are negligible. First we review the classical scattering theory 



in Sec. [II A| , before we discuss the quantum dynamics in Sec. [II B| . The longitudinal and 
transverse resistivities are discussed in Sec. [II B 3| , and resonance scattering is demonstrated 
in Sec. [II Q. The case of a 2DEG in a periodic array of fiux vortices, is the subject of Sec. |Il. 



In the first part. Sec. P^II A[ - pTD[ , of this section the general theory of electron motion in a 
periodic magnetic field is reviewed, and in the second part. Sec. [Ill E[ - [rTrF| , we present the 
results of the numerical calculations. 



11. SINGLE VORTEX SCATTERING 

In this section we will consider the introduction of magnetic flux vortices into a 2 dimen- 
sional electron gas, in the limit where each vortex can be treated as an individual scattering 
center. The vortices are assumed to be distributed at random, homogeneously over the sam- 
ple. The average separation between the vortices is assumed so large, that we can neglect 
interference from multiple scattering events. 

A. The Classical Cross Section 

We will start by calculating the differential cross section for a charged particle scattering 
on a flux vortex within the framework of classical mechanics. This will provide a reference 
frame, and allow us to speak unambiguously about the classical limit. In the calculations 
we will use an ideal vortex, which has a circular cross section with constant magnetic field 
inside, and zero magnetic field outside 

— 5- tor r < tq 
B{r) = I ^^0 (1) 

I for r > Tq. 

Here tq is the radius, and is the total flux carried by the vortex. The classical orbit is 
found as the solution to Newton's equation of motion with the force given by the Lorentz 
expression F = ev x B, for a particle of charge e. It consists, as is well known, of straight 
line segments outside the vortex, and an arc of a circle inside, with radius of curvature 
given by the cyclotron radius Ic = v/cuc, with v being the particle velocity, and cUc = — 
the cyclotron frequency. Because the orbit inside the vortex is an arc of a circle, and it is 
impossible to draw a circle that only cut the circumference of the vortex once, it follows 
that a particle which initially is outside the vortex, can never become trapped in the vortex. 
The geometry of the classical scattering process shown in Fig. |l], is determined by the 
dimensionless parameter 7 = Ic/tq, which is the ratio between the radius of the cyclotron 
orbit, and the radius of the flux vortex. The definitions of the reduced impact parameter 



P = b/ro, and the angles 0, ip, d, follows from Fig. |[ By inspecting the figure, it is observed 
that the following relations hold 
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where the sign of 7 is dictated by the direction of the magnetic field inside the vortex. Here- 
after we take 7 to be positive. After a small amount of arithmetic and if) are eliminated, 
and we have 



sm- = sign(7 + /3)W ,^,^^^, - 



(5) 



This relation gives the scattering angle 9 = 0{(3) as a function of the impact parameter jS. 
It is observed that for < 7 < 1 the scattering angle 6 sweeps through the interval [— tt, tt], 
for /3 G [-1,1], while for 7 > 1 we have 6 G [0, ^0], with Oq = 2arcsin(l/7) = ^(-1/7). 
Furthermore dd/d(3 = for /3 = — 1/7 producing a singularity in the differential scattering 
cross section at 6q. 

The differential cross section do jdQ gives the total weight of impact parameters, which 
give scattering into the direction Q. Equation |^ has at most two solutions, which are easily 
found to be 



/3^(^) = -7sin'0/2 ± cose/2Jl-72sin'^/2 



(6) 



Furthermore we have 
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The differential cross section is for < 7 < 1 given by 

dp. 



da 






de 

d(3+ 



dd 



for ^ < 
for ^ > 0, 



(8) 



and for 7 > 1, it is given by the expression 

da [--o^+ro^ tor 0<9<9, 



otiierwise. 

Examples of cross-sections and trajectories are shown in Fig. |^. The integrated cross section 
o"tot is equal to the total weight of impact parameters corresponding to particles which hit 
the vortex. It is equal to the diameter of the vortex citot = 2ro, as is always the case in 
classical scattering. If the classical cross section is substituted into the collision integral 
of the Boltzmann transport equation, it is straightforward as shown in App. 0, to obtain 
expressions for the transverse and longitudinal resistivities induced by a random distribution 
of vortices, in the limit where the density of vortices goes to zero, thereby eliminating the 
effects of multiple scattering events. Let B denote the average flux density applied to the 
sample, and p^f = B /n(.e the Hall resistivity of a 2DEG in a homogeneous magnetic field B. 
We then have 

Pxy _ /-clas Pxx ~ Pimp _ >clas /-irvX 

^xy 1 ^xx 1 \ ^) 

Ph Ph 

with 

The C,^^^'^ parameters are a function of the sole parameter 7 = Ic/tq. This is indeed what one 
would expect, since the transport coefficients calculated within a classical framework can 
not depend on the flux quantum 0o- The functional behavior of the (^-quantities, are shown 
in Fig. ^ We will postpone the discussion of the experimental consequences of these curves, 
until we have calculated the same quantities quantum mechanically. 

B. Quantum Dynamics 

In this section we shall consider the electron scattering on a magnetic flux vortex within 
the framework of quantum theory. We will calculate the electron wave function, from which 



the longitudinal and transverse conductivities can be found from the theory of Sec. [II B 3 . 
The quantum nature of the electron radically alters the picture of the scattering process, 
when the wavelength of the electron is comparable to, or longer than the diameter of the 
vortex. In the limit of very small electron wavelength, the scattering can essentially by 
described by the laws of geometrical optics, and thereby classical mechanics. 

1. The Scattering Wave Function 

We will again take an idealized cylindrical vortex, with constant magnetic field inside, 
and zero field outside, Eq. 0. This vortex is completely symmetric under any rotation about 
the center axis. This symmetry can also be made a symmetry of the Hamiltonian, by 
chosing the right gauge when writing down the vector potential. In cylindrical coordinates 
A = BrAr + egAg, we have 

B{r) = drAe - -dgAr + — . (13) 

When B is invariant under rotation, this equation has the simple solution 

Ar = 0, Ag{r) = - f dr'r'B{r'), (14) 

'^ Jo 

which in our case give the vector potential A = cgAg with 

{ — for r < ro 
^7 (15) 
for r > tq. 
27rr 

The Hamiltonian is given by the expression 

Here we have taken the charge of the particle to be e. The Hamiltonian is rotationally 
invariant, and therefore commutes with the angular momentum about the symmetry axis. 
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Lz- Consequently L^ and H have common eigenstates. The canonical momentum operator 
of a particle in a magnetic field is p = mv + qA = ^'V, and the operator for the angular 
momentum about the z-axis is L^ = [r x p]z = ^dg. The eigenstates of Lz are e*'^, where 
/ must be an integer in order not to have cuts in the wave function. We can now separate 
the variables of the common eigenstates of Lz and H, and write 

Mr,0) = RH{r)e''', (17) 

where k is an energy label E = f^k"^ /2m. Let us introduce the fiux quantum 0o = h/e 
and the dimensionless fraction a = 0/0o- The differential equations for the radial part of 
the wave function, takes a particularly simple form if we write them down in dimensionless 
variables ^ = t/tq and n = kr^. The variable k measures the size of the vortex compared 
to the electron wavelength. With these definitions, the equation for the radial part of the 
wave function for ,^ < 1 is 

i?" + -B! + (k" -{-- aif\ R = 0. (18) 



And for ^ > 1 we have 



^..(.-^^ 



R" + ±R'+{k'- ^^f^ )r = 0. (19) 



Inside the vortex an analytical solution to the radial equation can be found by following 
the procedure of L. Pagea, see also the review of Olariu and PopescuHiJ. The solution to the 
radial equation inside the vortex, which is regular as ^ — > 0, is 



')M 



{i,e) = ci\%~'^'^^'M 



2 



I la 



e'r (20) 



Here M is a confiuent hypergeometric function (solution to Kummer's equation), in the 
notation of Abramowits and Stegunli3, and C is a normalization constant which we will 
not need to evaluate. At the Landau quantization energies the confiuent hypergeometric 
function M reduces to an associated Laguerre polynomial. 

Outside the vortex the radial equation is just the differential equation for Bessel functions 
of the first kind. We have for ^ > 1 



Mi. 0) = (aiJi^^ikr) + biYi^^ikr)) e*'^. (21) 

The two constants ai, hi are found from the requirement, that the wave function has to be 
continuously differentiable at the boundary of the vortex. Let us write 

ai = ci cos 6i hi = -ci sin 6i, (22) 

then we have the following asymptotic expression for the total wave function ilj{r, 9) in the 
region kr '^ 1 

^{r, 6) = J-^ Yl C' co<^' + ^^y- (23) 

Here we have introduced the abbreviation xi = kr — (7r/2)(/ — a) — it /A. The scattering 
phase shift 6i, is found from the condition that the wave function has to be continuously 
differentially at the boundary of the vortex. The constant q can be expressed in terms 
of the phase shift, via the scattering condition, i.e. that all the incoming particle current 
must be described by the wave function ipi which far away from the vortex describes a 
uniform current directed towards the vortex. In the geometry with which we are concerned, 
the particles come in along the positive x-axis, towards the vortex placed at the origin, as 
shown in Fig. ^. The state describing a uniform current in the direction of the negative 
X-axis is 

^. = (,-ikr cos 9+iae_ (^24) 

This expression is observed to have a cut whenever a is not integer, however this cut will 
not be present in the total wave function ip. There will be a counter cut in the outgoing 
wave function ipg = ip — ipi, to render the total wave function completely smooth. This cut 
only arises because we insists in splitting the total wave function into incoming and outgoing 
terms. We take the angle 6 to have values in the interval — vr < 6 < tt, which places the 
cuts in ipi^s along the negative real axis. With this choice the asymptotic expansion of ipi, 
for /cr ^ 1, in terms of partial waves, takes the form 
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^„i„^ Z-^ L 



ikr—in/4 



cos7r(/ — a) + e 



'ikr+in/4 



]■ 



(25) 



V2iTkr ^ 

The scattering condition is now implemented by demanding that the coefficient of e"*'^'', in 
respectively the expansion oiip and ipi must be identical for each /, resulting in the identity 



ci = e 



i5i—in{l—a)/2 



(26) 



2. Phase Shifts for Scattering on Vortex with Finite Radius 



Let us now calculate the phase shifts for a charged particle scattering on a cylindrical 
flux vortex with finite radius. The phase shifts are derived from the condition that the 
logarithmic derivative of the radial wave function must be continuous at the boundary of 
the vortex 



1 dR< 



R< di 



1 dR> 



?=i 



R> di 



Inside the vortex we have with the abbreviations hi 



/| + l-/-(A;ro)V2«)/2and^/ 



^'-R< 



1 dRf 



di 



?=i 



hi M\hi + l,gi + l,a\ 

a + 2a -— -. . 

gi M[hi,gi,a\ 



(27) 

|/| + 1 
(28) 



Outside the vortex the logarithmic derivative reads 



1 dR> 



ji - yi tan 5i 



?=i 



Ji-Jk) -YjJK)tan6i' 



(29) 



where we have introduced the abbreviations zi = (A;ro/2)(Z/_Q_i(A:ro) — Zi_a+i{krQ)) with 
{z, Z) = (j, J) or {y, Y). It is now simple to solve for 6i 

ji - EiJi_aikro) 



ternSi 



(30) 



yi - EiYi^aikro)' 

The tan(5;'s are the bricks from which cross sections and transport coefficients can be build. 
The presented curves have all been calculated with phase shifts found by this expression with 
the help of Mathematical which have implementations of all the involved special functions. 
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3. Momentum Flow and Force Balance 

In order to express the transverse and longitudinal resistivities of the 2DEG in terms 
of matrix elements of the force, we make a simple force balance argument. For a more 
comprehensive discussion of the force balance, we refer to E. B. HansenEl. 

We consider a 2DEG carrying a current of density j^, along the x-direction, i.e. jy = 0. 
The distribution function of the electrons is then 

n{k) = n%k-6k), (31) 

with n°(fc) = 6(ei? — Cfc), and 6k is related to the current density jx = {evp/2h)g{eF)6k, 
where g is the density of states. The forces acting on an electron in the state k can be 
separated into two classes. First, there are the forces due to the electric fields E^, Ey 
present in the sample. Secondly, there are the forces due to scattering processes, which due 
to the rotational symmetry of the scatterers - vortices and impurities - can be decomposed 
into a force Fl parallel to fc, and a force Ft at right angles to fc. If we let 9^ denote the 
angle between k and the x-axis, we can write the scattering force on an electron in the state 
fc, projected along the coordinate axes as 

F,(fc) = FL{k) cos^fc - Erik) sin^^, (32) 

Ey{k) = EL{k) sin Ok + Erik) cos Ok- (33) 

When we have a steady state transport situation, the total force on the 2DEG must vanish. 
With n = x,y, and rie denoting the density of the 2DEG, we therefore have the equation 

ene^^ = 5^F^(fc)n(fc). (34) 

fc 

And thereby 

eUeEfj^ = -^g{tF)6k j —^E^{kF,9k)cos6k. (35) 

Along the x- and y-directins respectively, this leads to 
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^ = ^ (36) 

Ey Frikp) 



ja: ne^VF 



(37) 



The ordinary impurities etc. do not give rise to a net transverse force on the electrons. 
Consequently, only the vortices contribute to the transverse force ^^(fcir) on an electron, 
on the Fermi surface. This force will, in the limit of low vortex density, be proportional to 
the total number of vortices. In this limit we can therefore calculate the force as the total 
number of vortices N^-, times the force from a single vortex 



FT{kF) = N^{^>\Fy\^'j (38) 

Here Na = Aua = AB/a(j)o, with B equal to the average flux density applied to the sample. 
Further more, F = {e/2)[v x B + B x v] is the operator corresponding to the Lorentz force 
inside the vortex, and \1/ is the normalized wave function of the electron, i.e. (^I/|^) = 1. The 
electron wave function \1' is an eigenstate of the one electron - one vortex system. It describes 
an electron which initially travels along the x-axis, or in other words, the expectation value 
of the velocity is (\E'|'fi^|\l') = 6^xVf, and it is labeled by the k- vector k = {kp, 0). Note, that 
\1/ is identical to the wave function ip of Sec. |11 B 1| apart from the normalization {ijj\ip) = A. 
The transverse resistivity of the 2DEG can now be expressed as 

Pxy = — = PnCxy, (39) 

Jx 

where p'jj = B/riee is the Hall resistivity of electrons in a homogeneous magnetic field of 
magnetic flux density B, and 

Similarly, we have for the longitudinal resistivity 

Pxx - Pimp = PhCxx = ■:^-r'^^'^Cxx- (41) 



e^ n, 



e 



where we have singled out the contribution from the vortices, and 
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The matrix elements (tp\F^\-ip) can be calculated by numerically performing the integral, 
which is limited to the range of the magnetic field, i.e. to the core of the vortex. We have 
used this procedure, which is rather laborious for the computer, to check the results obtained 
by the following much faster procedure. 

In a continuum the flow of kinematical momentum is described by a continuity equation, 
which at the classical level takes the form 

dtk^ + d,7T,f, = F^. (43) 

Here k^ = mvn is the kinematical momentum (which in the presence of a magnetic field differ 
from the generalized momentum p^ = mf ^ + eA^) and tTj,^ = mVuV^ is the tensor describing 
the flow of kinematical momentum. The quantum mechanical equivalent of Eq. E3 is 



dtk, + d,7T,^ = ^*F^^, (44) 



where now 



k^ = m^*w^^, (45) 

171 

r^ut. = y ['^*{v.v^^) + {v^^){v,^r] , (46) 

and F is the operator corresponding to the Lorents force. The interpretation of Eq. ^ as 
a generalization of Eq. |4^ applies strictly speaking only to the real part. For a derivation 
of Eq. ^ we refer to the review of Olariu and Popescullll. The scattering situation we are 
considering is in a steady state and therefore dtk^ = 0. We have by Gauss law 

i/j\F^\^'j = Jd'rd,'K,,{r) = j rdd-Kr.ir), (47) 

where the index r refers to the radial component, and the last integral is along any circle of 
radius r > tq, (the Lorentz force is only non- vanishing inside the vortex). Let us first consider 
the y-component of the force. As 7r,y^ is a tensor, we can write TCry = iirr sin 6 + Tr^e cos 6. 
Furthermore the velocity operators are (outside the vortex core) 
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Vr = (48 

m I 

ve = —i--a), (49) 

mr I 



and due to the extra factor of r in ve, only 7r„ will contribute to the integral, Eq. ^ in the 
limit r — i> oo. In the asymptotic region we have 



TTr 



VTTTir -^ — ^ -^ — ^ Z 



^ ^ ei(5,+„-50-if m cos[5i+„ - 5/ - -mY'^\ (50) 

and 

>|FJ^\ = / rd^TT^^sin^, (51) 



resulting after some algebra in 

^-J' = ^E^i^[2(5,+i-50]. (52) 

Similarly, we have in the x-direction -Krx = ^Crr cos 6 — lire sin 6*, and 

U\F^\iIj\ = I rdeTfrrCose, (53) 

J — 7r 

resulting in the expression 

Cxx = — V sin^ [5^+1 - 6i] . (54) 

Tia ^ — ^ 
I 

In order to calculate the ^-quantities it is convenient to express them in terms of ti = tan Si 



{ti+i-ti){i + titi+i) 



^ 1 ^ {ti+i-ti){l + tiU 



(56) 



The Ca;.j^ expression can be divided into two terms of different origin. Define ^Cxy and ^xy{o) 
to be respectively 

2 ^^ {U+i — U)titi 
I 

A (^\-r A/- (t;+i -ti){l -titi+i) , . 

A.,(«) = C.y - AC, - (l + t2)(l+t2_^J • (^^^ 
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Then we have 






l + tl, 1 + tf 



= — [sin 600 cos 600 — sin 6-od cos 6-00] 
Tia 

= '^^^^. (59) 

Here we have used the fact that the sum is telescopic, and the following properties of the 
phase shifts. For any given finite energy of the electron, the Tth partial wave, which is an 
eigenstate of L^ with eigenvalue hi, will have an insignificantly small propability density 
inside a radius 6/, determined by h\l\ ~ hihk, i.e. hi ~ \l\/k. As |/| grows to infinity the Z'th 
partial wave will not be able to detect any difference between a vortex of radius vq <^ bi, and 
an Aharonov-Bohm string, and consequently Si -^ 5f^ for / -^ ±00. The Aharonov-Bohm 
phase shifts are 

I for / > a 
6^ = { (60) 

7r(/ — a) for / < a, 



which can be seen by comparing the wave functions given respectively in Eq. ^ and Eq. pl| . 
By inserting the Aharonov-Bohm phase shift into Eq. ^, the result follows 

_ Sin27ra (t;+i - t;)t;t;+i 

By inserting 6f^^ into A(.j.y, it is seen that Ai^^y vanishes in this limit. In this form (^y is 
therefore separated into a term due to the Aharonov-Bohm effect, and a correction term due 
to the effects of the finite radius of the vortex. Moreover the sum over / in the correction term, 
converges rapidly for increasing \l\. The A^y{a) term show, that even though the magnetic 
field is only non- vanishing at a single point in the Aharonov-Bohm limit, the electron still 
experiences a transverse force. This force was first described by S. V. lordanskiiEJ. Let us 
remark that the sum appearing in (^x is already convergent. The residual value of (xx in the 
Aharonov-Bohm limit is 

e/ = i-sin2vra. (62) 
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Curves showing Cxy as a function of kr^ for different values of the flux a, have been plotted 
in Fig. ^. For kr^ ^ 1 all the curves go to C,xy = 1, which is the classical value as is seen in 
Fig. 1^. In the Aharonov-Bohm limit kr^ <^ 1, the curves go to a residual value depending 
on the flux (^y -^ Cx^ = sm27ra/27Ta. The Cxy-curve for a = 1/2 we can compare with the 
experimental Hall factor measured by Geim et al.a. The overall qualitative behavior is in 
good agreement. To make a quantitative test we have fitted our curve to the experimental 
curve by tuning the radius of the vortex tq. The best fit is obtained for a vortex radius 
tq = 50nm, and this is significantly smaller than the exponential length estimated by Geim 
to be lOOnm. 

Curves showing 27ra(xx as a function of kvQ for different values of the flux a, is shown in 
Fig. Ij. The magneto-resistivity, given by 

^ = -C.., (63) 

dB ne 

is constant, corresponding to the fact that the resistivity is proportional to the density of 
vortices, in the low 5-field limit. The magnitude of the magneto-resistivity do not show 
perfect agreement with the values measured by Geim et al.S. For the 2DEG with the lowest 
density that has been measured by Geim, ng = 3.65 ■ lO^'^cm"^, Geim reports dp^^^ /dB = 
O.OQQ/G. In this case Eq. ^ predicts dpxx/dB = 0.55Q/G, for tq = 50nm. 

In order to investigate the sensitivity of the C's to the radial distribution of magnetic 
field in the vortex, we have calculated them with an exponential distribution 

Be.Ar) = l^e-^/"- (64) 

The calculation is performed by solving the radial Shrodinger equation numerically, and im- 
posing an artificial boundary on the vortex, outside which the amount of flux is negligible. 
The phase shifts are then found from the condition that the wave function much be smooth 
at the boundary. From the phase shifts the ^-quantities can be calculated as for the step 
function vortex. The calculation shows almost no difference between the ^-values, corre- 
sponding to the two different choices, indicating that the ^-quantities are not very sensitive 
to the shape of the vortices. 
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C. Multi Flux Quantum Vortex and Resonance Scattering 

When the total amount of magnetic flux inside the vortex is increased, the ^-quantities 
acquire more structure. To illustrate this we have plotted (xy and 27ia(xx in Fig. for a flux 
vortex carrying a total of a =10 flux quanta. The structure seen in the plots is an effect of 
the resonant scattering which takes place when the energy of the incoming particle is close to 
one of the Landau quantization energies corresponding to the magnetic fleld strength inside 
the vortex. The magnetic fleld in the vortex vanishes outside a flnite range - the radius of the 
vortex - and there are therefore no real Landau levels in the sense of stationary eigenstates, 
but only metastable states. The resonance condition {p + l/2)tujjc = h^k'^/2m translates for 
flxed a into 



Vo = 2y|a|(p+^), p = 0,l,2,... (65) 

These values are in excellent agreement with the resonances seen in Fig. |^, where the flrst 
eight resonances corresponding to p = 0, . . . , 7, are clearly distinguished. At the resonance 
energies the typical time the particle spends in the scattering region, i.e. inside the vortex, is 
much longer than it is away from resonance. The time the particle spends inside the vortex 
at a resonance, is the lifetime of the corresponding metastable state. The inverse lifetime is 
proportional to the width of the resonance, that is strictly speaking the width of the peak in 
the partial wave cross section ai, corresponding to the / quantum number of the metastable 
Landau state. 

It is easy to interpret the peaks in the Cxx-curve in Fig. |^, appearing at the resonance 
energies. Because when the electron spends longer time in the scattering region, it loses 
knowledge of where it came from, resulting in an enhanced probability of being scattered in 
the backwards direction. The Cxy-curve is a measure of asymmetric scattering, and we can 
therefore interpret the dips seen in Fig. |^ along the same line of reasoning. The electron 
spends longer time in the scattering region, thereby losing knowledge of what is left and 
what is right. Every time a new "channel" is opened the asymmetry of the scattering is 
suppressed, thus resulting in a sawtooth like curve. 



III. HALL EFFECT IN A REGULAR ARRAY OF FLUX VORTICES 

A. Introduction 

Recently measurements were made by Geim et al.cfa of the Hall resistivity of low density 
2DEG's in a random distribution of flux vortices, at very low magnetic field strengths. A 
profound suppression of the Hall resistivity was found, for 2DEG's with Fermi wavelengths 
of the same order of magnitude as the diameter of the flux vortices. This indicates that 
we are dealing with a phenomenon of quantum nature. These measurements were made by 
placing a thin lead film on top of a GaAs/GaAlAs heterostructure. When a perpendicular 
magnetic field is applied, the magnetic field penetrates the superconducting lead film and 
also the heterostructure in the form of flux vortices each carrying half a flux quantum 0o/2 
of magnetic flux. Due to the strong flux vortex pinning in the films Geim have used, the 
vortices were positioned in a random configuration. 

In this section we will consider the hypothetical experiment where one instead of a "dirty" 
film, places a perfectly homogeneous type II superconducting film, on top of the 2DEG. The 
film do not have to be made of a material which is type II superconducting in bulk form. 
A film of a type I superconducting material will also display a mixed state if the thickness 
of the film is below the critical thickness dc- Experimentally perfect Abrikosov fiux vortex 
lattices have been observed in thin films of lead with thickness d < dc = 0.1/im,t3. 

If one succeeds to make such a sandwich construction, one has an ideal system for 
investigating how a 2 dimensional electron gas behaves in a periodic magnetic field. When 
the magnetic field exceeds Hd, which can be extremely low, the superconducting film will 
enter the mixed phase, and form an Abrikosov lattice of fiux vortices. The Abrikosov lattice 
in the superconductor will give rise to a periodic magnetic field at the 2DEG, and moreover 
as the strength of the applied magnetic field is varied the only difference at the 2DEG, is that 
the lattice constant of the periodic magnetic field varies. The Abrikosov lattice is most often 
a triangular lattice with hexagonal symmetry, although other lattices have been observed 
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(e.g. square) in special cases where the atomic lattice structure impose a symmetry on the 
flux lattice,Ej. For simplicity, the model calculations we have made are for a system with 
a square lattice of flux vortices, but we do not expect this to influence the overall features 
of the results. From the point of view of the 2DEG it is important that the flux vortices 
carry half a flux quantum h/2e = 0o/2 due to the 2e-charge of the Cooper pairs in the 
superconductor. The magnetic field from a single flux vortex fall of exponentially with the 
distance from the center of the vortex. This exponential decay is characterized by a length 
As, which essentially is the London length of the superconductor, proportional to one over 
the square root of the density of Cooper pairs. The length As can be varied by changing the 
temperature, or the material of the superconductor. 

The other characteristic lengths of the system are the Fermi wavelength Xp = A/27r/r7^, 
where rig is the density of the 2DEG, the lattice constant a of the periodic magnetic field, 
and the mean free path If = vpT. The mean free path we assume to be very large compared 
to a and As. The magnetic field is only varying appreciably when a is larger than As. This 
means that the magnetic flux density of the applied field should be appreciably less than 
0o/(7rA^). In the limit where As <^ a,Xp, the vortices can be considered magnetic strings, 
and the electrons experiences a periodic array of Aharonov-Bohm scatterers. In this case 
the value of the flux through each vortex is crucial. If for instance the flux had been one flux 
quantum 0o, the electrons would not have been able to feel the vortices at all. But in the 
real world the vortices from the superconductor carry 0o/2 of flux, and therefore this limit 
is nontrivial. The electrons has for instance a band structure quite different from that of 
free electrons. In the mathematical limit of infinitely thin vortices each carrying half a flux 
quantum, there cannot be any Hall effect. This is most easily seen by subtracting one flux 
quantum from each vortex to obtain a flux equal to minus a half flux quantum through each 
vortex. As we have discussed earlier the introduction of the Dirac strings can not change 
any physics, and the procedure therefore shows that the system is equivalent to it's time 
reversed counterpart, thereby eliminating the possibility of a Hall effect. 

In this study we have ignored electron-electron interaction and the electron spin through- 
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out, in order to keep the model simple. From the point of view of the phenomena we are 
going to describe, the effect of the electron spin will be to add various small corrections. 

B. Electrons in a periodic magnetic field 

1. Magnetic translations 

It is a general result for a charged particle in a spatially periodic magnetic field B{x,y), 
that the eigenstates of the system can be labeled by Bloch vectors taken from a Brillouin 
zone, if and only if the flux through the unit cell of the magnetic fleld is a rational number 
p/q times the flux quantum. The standard argument for this fact is made by introducing 
magnetic translation operators. To introduce magnetic translation operators in an inho- 
mogeneous magnetic fleld we flrst make the following observation. The periodicity of the 
fleld can be stated B{r -\- R) = B{r), for R belonging to a Bravais lattice. But this im- 
plies that the difference between the vector potentials A{r + R) and A{r) must be a gauge 
transformation 

V X {A{r + R)- A{r)} = B{r + R) - B{r) = 0. (66) 

We introduce the gauge potential xr ^^^ write 

A{r + R) = A{r) + VxRir). (67) 

The function xr is only deflned modulo an arbitrary additive constant which have no phys- 
ical effect. The Hamiltonian of the electrons is 

H = ^{p-eA{r))\ (68) 

The ordinary translation operators Tr = exp[j^R-p] do not commute with the Hamiltonian, 
because they shift the argument of the vector potential from r to r + R, but as we just have 
seen this can be undone with a gauge transformation. We therefore introduce the magnetic 
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translation operators as the combined symmetry operation of an ordinary translation and a 
gauge transformation 

6 I 

Mfi = exp[-2-Xii(r)]exp[-i2-p]. (69) 

The operator Mji is unitary, as it is the product of two unitary operators, and therefore has 
eigenvalues of the form e*'^. Let us denote the primitive vectors of the Bravais lattice a and 
h. We can find common eigenstates of Ma-, Mb and H, if and only if they all commute with 
each other. The magnetic translations each commute with the Hamiltonian by construction, 
and furthermore we have 

MaMb = exp [2m^] M^M^. (70) 

If the fiux through the unit cell is a rational number p/q (p and q relatively prime) times 
the fiux quantum (pQ, Mqa and M^ commute. In this case the cell spanned by qa and h 
is called the magnetic unit cell. Let us define c = qa. The possible eigenvalues of Mf. are 
phases e^'^*'^i, where we can restrict |A;i| < |, and equivalently for M5. We can therefore label 
the common eigenstates |A;, n), where k = kic* + k2b*, and c*, b* are the primitive vectors of 
the reciprocal lattice. The vector k is restricted to the magnetic Brillouin zone. An arbitrary 
magnetic translation of an eigenstate with a Bravais lattice vector R = nc + mb can now be 
written Mji\k,n) = (Mc)"(M5)™'|fc,n) = exp[ik ■ R]\k,n), showing that the eigenstate is a 
Bloch state. In this case we can speak of energy bands forming a band structure in the usual 
sense. When the fiux through the elementary unit cell {a, b) is an irrational number of fiux 
quanta, the situation is different. The irrational number can be reached as the limit where 
p and q get very large, and consequently the Brillouin zone get very small and collapses in 
the limit. 

2. The Dirac vortex viewpoint 

In this section we will show, how it is possible to argue in a slightly different way from 
the previous section, and hereby in a simpler way obtain the vector potential of a periodic 
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magnetic field. Let us again assume a rectangular unit cell (a, 6), B{x + a,y) = B{x,y + 
b) = B{x,y) etc., to keep the notation simple. The magnetic field enter the Hamiltonian 
only through the vector potential. The question is therefore if one can choose a gauge 
such that the vector potential will be translationally invariant relative to a unit cell (c, d) 
A{x + c, y) = A{x,y + d) = A{x, y) etc. It is clear that if such a periodic A-field exists, then 
the total fiux $cd through the unit cell (c, d) will be zero, as it is given by the line integral 
of A around the boundary of the unit cell (c, d) 



^cd= j> A dl, (71) 

Jd{c,d) 

which is zero by the periodicity. We remark that due to the relation B = 'V x A, the cell 
(c, d) will be bigger than or equal to (a, b). If the fiux through the unit cell of the magnetic 
field is not zero, but equal to a rational number times the fiux quantum, ^ab = p/q <Po ip 
and q relatively prime), a trick can be applied to make the fiux $cd become zero. It is a 
basic fact, apparently first observed by DiractH, that a particle with charge e cannot feel 
an infinitely thin solenoid carrying a fiux equal to an integer multiple of the fiux quantum 
00 = h/e. Such a stringlike object carrying one fiux quantum is sometimes called a Dirac 
vortex. To find the periodic vector potential take an enlarged unit cell (c, rf) = {qa,b), so 
that $cd = P0O and put by hand p counter Dirac vortices through the cell, to obtain zero net 
fiux. Then a divergence free vector potential can be build for instance by Fourier transform 

B{Q) = 4 / rfVexp[-iQ ■ r]B{r), (72) 

cd J(c,d) 



(c,d) 

\ 



A{r) = j:' '"^^ 
Q^o \ -iQ, 



J 



'-^^^^BiQl (73) 



where the sum is over Q in the reciprocal lattice. Here we have used continuum notation, 
but it is straightforward to write down the lattice equivalents of the expressions. 
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C. Lattice calculation of Hall conductivity 

We have calculated the Hall conductance of the 2DEG in the vortex field by a numerical 
lattice method. This we do because the calculations then reduces to linear algebra operations 
on finite size matrices, which can be implemented in a program on a computer. The idea 
is to consider an electron moving on a discrete lattice, rather than in continuum space. 
We know, although we are not going to prove it here, that in the limit where the discrete 
lattice becomes fine-grained compared to all other characteristic length of the system, the 
continuum theory is recovered. Here we assume that the original Bravais lattice has square 
lattice symmetry, with a lattice parameter which we call a. The discrete micro lattice is 
then introduced as a fine-grained square lattice inside the unit cell of the Bravais lattice. 
The lattice parameter of the micro lattice we take as a/rf, where d is some large number, in 
order to keep the two lattices commensurable. The condition that we have to impose on the 
micro lattice, in order that it is a good approximation to the continuum, can then be stated 

a/d <ti a, Xp-, As, • • • . (74) 

In the numerical calculations we have made, we have taken d = 10. 
The tight-binding calculations are made with the Hamiltonian 

H = -Yl ^i+^.J+^'4+xS-+x'- (75) 

iJTT-' 

Here i, j are Bravais lattice vectors, and r, r' are vectors indicating the sites in the basis. 
The matrix elements ti^Tj^T' are taken non-zero only between nearest neighbour sites. The 
matrix element between two nearest neighbour sites r and r' are complex variables t-r^T+Bf,, = 
j-^tA^{T) ^]^\^ g^ phase given by the vector potential A^{t) residing on the link joining the sites. 
The translation invariance of the Hamiltonian can then be stated ti+T^ij^T'^i = ti^Tj+T', 
for all vectors I belonging to the Bravais lattice. Let us introduce the system on which our 
calculations were made as an example. Fig. ^ shows the unit cell with its internal structure 
i.e. the basis. There are N = d ■ d sites in the basis. The length of the links we write 
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as a/d, where a is the side of the unit cell, with area Vt = a^. The vectors r = (ri,r2)^, 
Ti,T2 = 0, l,...(i — 1 are offsets into the basis, while the vectors i = {ii,i2)ci, ii,i2 G ^ 
indicate the cells in the Bravais lattice. The operator ct^^, for a given r, is defined on the 
Bravais lattice, and accordingly it can be resolved as a Fourier integral over the Brillouin 
zone as 



ct 



Q 



rf2^g-iq.y+x) +^_ (76) 



Inserting this and using the translation invariance, the Hamiltonian can be rewritten as 

where we have introduced 

Hk = -J2 ^xJ+x'e'''^^'+"'""^4,xCfc,x'. (78) 

JT-r' 

It is seen that Hk only mixes the A^ states jfer), i.e. it is a A^ x A^ matrix. The A^ eigenvalues 
of i/fc are the energies of the A^ tight binding Bloch states with wave vector k. Let us denote 
the eigenstates of Hk by -u^ 

H^ul = Elul (79) 

where a = 1, 2, ... A^ and iJ^" < E'^'^^. From the A^ dimensional vector -u^ we can construct 
the eigenstate '^% of the Hamiltonian H 

(j + r|^^) = e*(^'+")<(r). (80) 

It is straightforward to verify that this is the correct Bloch eigenstate of H . The band struc- 
ture can be calculated directly by diagonalizing the N x N matrices, H^, for representative 
choices of k in the Brillouin zone. Before one can compare the spectrum obtained from 
this calculation with that of a continuum system, a scaling of the energies is required. To 
scale the energy to the spectrum of a particle with an effective mass m, we have to take 
t = h^(P/ma^, and e„(fc) = E'^ + At. 
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The Hall conductivity can be calculated by the same method as in the homogeneous 
magnetic fieldQ'Ea. We have used a single particle Kubo formula to calculate the Hall con- 
ductance 

— i^ \^ {Jx)ali{Jy)l3a - {Jy)ap{Jx)pa ,^.s 

E°'<Ef<EP ^ ' 

where Jx-,Jy are the currents in the x^y directions, and the sum is over single particle states 
|a, h) with energies below and above the Fermi level Ep. The area of the system is denoted 
Aq. All quantities are diagonal in fc, and therefore this index is suppressed. The summation 
is composed of a discrete sum over bands, and an integral over the Brillouin zone for each 
band. The Brillouin zone shown in Fig. ^ is doubly connected because the states on the 
edges is to be identified according to the translation invariance. This gives the Brillouin zone 
the topology of a torus T^, with two basic non-contractible loops. The current operator can 
be written 

^ ' rf'fcjfc, (82) 

BZ 



(27r)2 



where J ^ = ie/h){dHk/dk). By use of some simple manipulations and completeness, it is 
straightforward to refrase Eq. |8l| 



ie sr^ I / da 



E'^<Ep 



"^ hAo „^„ V \ dt,.. 



da \ I da 



C/C-y / \ UfX/y 



da 



dh: 



(83) 



where 



— ^ is shorthand for ^|/5, fc). This formula was first derived by Thouless, 



dkfj, I dk^ 

Kohmoto, Nightingale and den NijsQ, for a noninteracting two-dimensional electron gas 
in a periodic scalar potential, and a commensurate perpendicular magnetic field. It requires 

M^ ) it is necessary to consider the 



some comments to be meaningful. In order to calculate 
difference {\a,k + 6k^) — \a, k))/6k^. But this difference is not well defined as it stands, as 
the phase of the states is arbitrary. Rather than representing the state Uk by a single vector 
in C^, it should be represented by a class of vectors which differ one from another only by 
a phase. These equivalence classes are sometimes called rays. To compare states locally, we 
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need to project this U{1) degree of freedom out. This is done by demanding the wave func- 
tion to be real and positive, when evaluated in a fixed point, i.e. u'^{Ti) = {Ti\a, k) G M_|_. 
If the wave function happens to be zero in Tj, some other point Tj must be used. When a 
band has a non-zero Hall conductivity, it is not possible to find a single r which work for 
all the states in the Brillouin zone. The change from Tj to Tj which shifts the phase of the 
states, is analogous to a "gauge transformation" on the Brillouin zone, of the set of states. 
The special combination of terms which appear in Eq. ^ is gauge invariant with respect 
to these special "gauge transformations". If we let |x") denote a state which is obtained 
from I a, k) by fixing the phase according to the above scheme, the following formula for the 
contribution to the Hall conductivity from a single band a, is well defined 



,;„ = £! 1/ ,=J/«X° 



^•^ h 2m Jbz I \ dky 



dx'^ \ I dy 



dx" 



It has been shown in detail by Kohmotoc3, for the homogeneous magnetic field case, that 
this expression is equal to minus e^ jh times the first Chern number of a principal fiber 
bundle over the torus. As the first Chern number is always an integer, this has the physical 
consequence that whenever the Fermi energy lies in an energy gap, the Hall conductance 
is quantized. We will use this result to interpret certain peaks in the o"a;y-spectra we have 
calculated. When the Fermi level is not in an energy gap of the system, we will have to use 
Eq. |8^ to calculate a^y. As we shall see, in this case there is no topological quantization of 
the Hall conductivity. 

D. Energy band crossing 

In this section we study the effect on the Hall conductivity of an energy band crossing. 
This has previously been discussed in different contexts by several authorsi^Ej. 

When the shape of the magnetic field is varied, controlled by some outer parameter ^, 
it will happen for certain parameter values ^o; that two bands cross, see Fig. 110. This is 



the consequence of the Wigner-von Neumann theorem, which states that three parameters 
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are required in the Hamiltonian in order to produce a degeneracy not related to symmetry. 
Here the parameters are kx,ky and the outer parameter ^, which in our calculation is the 
exponential length of the flux vortices from the superconductor. When the energy difference 
E~^ — E~ between the two bands considered is much smaller than the energy distance to the 
other bands, the Hamiltonian can be restricted to the subspace spanned by the two states 
|+, fc") and |— , kP). The point in the Brillouin zone where the degeneracy occur we denote 
fc°. The Hamiltomian H{k) is diagonal for k = fc°, and we denote the diagonal elements 
respectively Eo + e and E^ — e. For small deviations of k from fc" the lowest order corrections 
to the Hamiltonian is off-diagonal elements A(fc) linear in A; — k^ . Without essential loss 
of generality we can assume that e is independent of k. Then e plays the role of the outer 
parameter controlling the band crossing. The Hamiltonian is then approximated by 

/ 



H{k) 



\ 




(85) 



The off-diagonal element is expanded as 

A{k)=a{k^-kl)+P{ky-kl) (86) 

where we have introduced the matrix elements a = ^(— , fc°|iJfc|+, fc°), and P = 

We want to find the consequences of the energy band degeneracy, on the topological Hall 
quantum numbers of the bands. Let us define 



^-C^' H ( |± 



d±\ / d± 



m 



dky 



Then the interesting quantities are the integrals of B^{k) and B^{k), around a small neigh- 
bourhood of the degeneracy point k^. It turns out, that it is the two numbers a and f3 that 
control what happens. 

In general a and (3 will be nonzero complex numbers — nonzero because we have assumed 
the degeneracy to be of first order. Let us first consider the degenerate case where a and (3 
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are linearly dependent, i.e. a/ (3 is real, or otherwise stated lm.{a*l3) = 0. Then by a linear 
transformation we can write the Hamiltonian 



H(k) 



(ki + K,2)a^ + ea^ 



e Ki + K2 



V 



Kl + /«2 



(88) 



/ 



where Ki, k^ are rescaled momentum variables, with dimension of energy. The cr^'s refers to 
the Pauli matrices. The eigenvectors of this Hamiltonian can be chosen real, and this will 
clearly make BjJ^k) vanish. In this case we therefore conclude that there is no exchange of 
topological charge. Here the word topological charge is used to denote Hall quanta. 

Let us now treat the general case where a and /9 are linearly independent, i.e. Im(Q;*/3) ^ 
0. In this case we can, by a linear transformation, write A(fc) = /ti + ik^ = /te*^, which 
defines the scaled variables k, 6. This reduces the Hamiltonian to the form 



/ 



H{k, 9) = Kia^ + K20-^ + ea^ 



KC 



\ 









(89) 



Let us define A = veM-fi?. Then the eigenvalues of H{k,, 6) are ±A and the two corre- 
sponding eigenstates are 



/ 



|±,fc) 



v^ 



±^l±e/\ 



\ 



\ 



e^Vl^e/A 



(90) 



In order to calculate the integrals of the i?±-functions we need to express them in terms of 
the K, ^-variables. The Jacobian of the transformation is given by the expression 

ndndO 



and 



B±{k, 



UjKxOj'^11 



CX'pPi CxiP'p 



\(lfPi CXiP'p 



(91) 



d± 



d± 



d± 

dn 



d± 



(92) 



where the indices r, i refer to the real and imaginary parts respectively. We can now calculate 
the contribution to the Hall conductivity from each of the bands, from the area around k'^ 
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given by |fi;| < Kc, where Kc is some local cutoff parameter which limit the integration to the 
area where the approximation leading to the Hamiltonian Eq. ^ is valid 

,2 



:sign[Im(a*/3)] / dK dO 
' Jo Jo 



h2m ^ ' ' "^'^ Jq Jo 2A3 

^ . . .„, e Z"""^'"' udu 



±^sign[a*/5]-y^ (1 



+ 7/2)3/2 



= ±^sign[Im(a*/5)]sign[e], (93) 

where the last equality sign is valid when Hc/\^\ ^ 1, i-e. close to the crossing where e = 
0. Here the factor sign[e] signals that the two bands exchange exactly one topological 
conductivity quantum e^/h, at the crossing. This is not surprising, because we know that 
the total contribution from the states in a single band is always an integer times the quantum 
e^/h. Moreover it is readily seen that in the hypothetical situation of a n'th order degeneracy, 
i.e. one for which A = /t^e*"^, n quanta are exchanged. The total topological charge in the 
band structure is conserved. The topological charge can flow around and rearrange itself 
inside a band, but only be exchanged between bands in lumps equal to an integer multiplum 
of the conductivity quantum e^/h. When we gradually shrink the radius of the flux vortices 
to zero, the Hall effect has to disappear. There are two mechanisms with which the Hall 
effect can be eliminated. The first is by moving the topological charge up through the band 
structure by exchanging quanta, resulting in a net upward current of topological charge, 
eventually moving the charge up above the Fermi surface, where it has no effect. The 
second mechanism is by rearranging the topological charge inside the bands, so that each 
band has a large negative charge in the bottom, and a large positive charge in the top, 
but arranged in such a clever way that charge neutrality is more or less retained for all 
energies. This second mechanism will also give a net displacement of topological charge up 
above the Fermi energy, because in general the Fermi surface cuts a great many bands, and 
for all these bands the large negative charge, which they have in their bottom part, will 
be uncompensated. To use the language of electricity theory, we can say that every band 
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gets extremely polarized, resulting in a net upward displacement current, in analogy with 
the situation in a strongly polarized dielectric. Our numerical calculations indicate, that it 
is the second mechanism which is responsible for the elimination of the Hall effect, as the 
radius of the vortices shrinks to zero. 

When two bands are nearly degenerate for some fc°, each of the bands have concentrated 
half a quantum in a small area in k-space around fc°, and in general the topological charge 
piles up across local and global gaps in the energy spectrum. This is the reason for the 
oscillatory and spiky behavior of the Hall conductance as a function of electron density, that 
is seen on the calculated spectra below. It is also the reason why the numerical integration 
involved in the actual evaluation of the Hall conductivity, is more tricky than one could 
wish. 

E. Numerical results 

1. Transverse Conductivity 

We have calculated the transverse conductivity axy as a function of the integrated density 
of states, for electrons in a square lattice of flux vortices, for a series of varying cross sectional 
shapes of the flux vortices. Each of the field configurations consists of a square lattice of 
flux vortices with a given exponential length A^. The parameter which vary from calculation 
to calculation, is the dimensionless ratio ^ = A^/a, where a is the length of the edges of the 
quadratic unit cell. The unit cell is shown on Fig. ^jand contain, as we have already discussed, 
two vortices and a counter Dirac vortex. In order to do the tight binding calculation a micro 
lattice is introduced in the unit cell. In all the numerical calculations we present, the micro 
lattice is 10 x 10. This gives 100 energy bands distributed symmetrically about the center 
on an energy scale. Out of these only the lower part, say band 1 to 20, approximate the real 
energy bands well, while the rest is significantly affected by the finite size of the micro lattice. 
A careful examination of the vector potential reveal that the symmetry of the Hamiltonian 
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is very high for the particular choice of unit cell shown in Fig. ^. The field from a single 
vortex we have taken as i^Qe^^l'^^^l+l'^^'l)/''''' instead of the more realistic Bqc^^'^^^'^^ With this 
choice the vector potential can be written down analytically in closed form. This makes the 
calculations simpler, and does not break any symmetry that is not already broken by the 
introduction of the micro lattice. (We have made calculations of the band structure, with 
both kinds of flux vortices, and the differences are indeed very small). The energy spectrum 
is invariant under the changes {k^, ky) i — > (iA;^, ^ky), {^ky, ifc^^). This fact is exploited to 
present the band structures in an economic way. The labels F, X and M correspond to the 
indicated points in the Brillouin zone. Fig. p. 



In Fig. [TT] we have plotted a selection of typical band structures which illustrates the 
crossover from the completely fiat Landau bands in the homogeneous magnetic field ^ = 00, 
to the band structure of electrons in a square lattice of Aharonov-Bohm scatterers with 
a = 1/2 at ^ = 0. The band structures have been found by direct numerical diagonalization 



of the Hamiltonian. (See also Fig. \l^ . 



In Fig. O some typical results of the numerical calculations of the Hall conductivity 



are shown. In general we have no reason to expect, that the Hall conductivity should be 
isotropic as a function of the angle between the current and the flux vortex lattice. The 
results we present is for a current running along the diagonal of the square lattice, i.e. along 
the X-axis in Fig. |^. It is seen, that whenever there is a gap in the spectrum, the Hall 
conductivity gives the quantized value in agreement with the discussion in the last section. 
At Fermi energies not lying in a gap a^y always tend to be lower than the value it has in 
the homogeneous field. And in the limit ^ 1 — > 0, a^y vanishes altogether. In this limit 
the electron sees a periodic array of Aharonov-Bohm scatterers, each carrying half a flux 
quantum, and there is no preferential scattering to either side. We observe that when the 
flat Landau bands starts to get dispersion, the contribution to the Hall effect is no longer 
distributed equally in the Brillouin zone. Instead it piles up across local and global gaps in 
the spectrum resulting in the spiky a^y spectra. The density of topological charge is plotted 



as a function of the fllling fraction in Fig. |13| It is seen that for ^ ^ the distribution gets 
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strongly polarized, with a negative contribution to the Hall effect at the bottom part of the 
bands, and a positive contribution at the uppermost part of the bands. 

In Fig. ^, the filling fraction is limited to values below 10, but the same behavior 
continues for larger values of the filling fraction. We have observed the spiky behavior up 
to z/ = 30, and have no reason to believe that it should not continue to larger values. The 
limitations on our calculations, comes from the fact, that we are only able to handle matrices 
of limited size in the numerical calculations. Filling fractions of order 30 correspond in this 
context, to very low density electron gases, and consequently our calculations belong to the 
"quantum" regime, i.e. to the regime where Xp ^ A^. According to the discussion in Sec. ||, 
we expect a cross-over to a semiclassical regime, for electron gases with higher density, where 
Xp '^ Xs, with a qualitatively different behavior. 

It is an important question whether it is possible to observe these features of the Hall 
conductivity in experiments. The conditions, which are necessary, are that the mean free 
path I is long compared to all other lengths, and that a, Xg < Xp. The last condition 
indicates that we are concerned with quantum magneto transport, in the sense that the 
vector potential is included in the proper quantum treatment of the electrons. This is 
in contrast to the case Xp <^ a, Xg where the electron transport can be treated as the 
semiclassical motion of localized wave packets in a slowly varying magnetic field. Let us 
assume, in order to make some estimates, that the superconductor has a London length 
somewhat less than lOOnm, resulting in an exponential length of the vortices A^ of about 
lOOnm at the 2DEG, after the broadening due to the distance between the superconductor 
and the 2DEG has been taken into account. In order to have a variation in the magnetic 
field we should have a > Xg, and to be in the quantum regime A^ > a, Xg. This gives the 
estimate for the electron density Uf. ~ 10^°cm~^, which is not unrealistic. The effect of 
the impurities is (to first order), to give the electrons a finite lifetime. This gives a finite 
longitudinal conductivity axx and broadens the density of states. It also introduces localized 
states at the band edges (Lifshitz tails). If the field is nearly homogeneous, we have the 
standard quantum Hall picture with mobility edges above and below every Landau band, 
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resulting in the formation of plateaus in a^y, which only can be observed at much higher 
magnetic fields, of order lO^G, where the filling fraction is of order one. On the other hand 
when the amount of impurities is low, that is kpl ^ 1, all these effects will be small, and 



we expect that the features of axy, shown in Fig. 12 will have observable consequences in 



Pxy = (^xyKcrlx + (^"^ 



xyJ 



In Fig. |1J we have plotted the band structure of electrons in a square lattice of vortices 



carrying one flux quantum each. The calculation has been made with the same basis as the 



band structures shown in Fig. |ll], the only difference is that all fluxes have been multiplied 
by a factor of 2, as can be seen from the spacing between the Landau bands. In the limit 
where ^ — >■ 0, and the vortex lattice becomes a regular array of Aharonov-Bohm scatterers, 
we recover the well known band structure of free electrons. This is in agreement with our 
discussion of the AB- vortex in Sec. ||. 

2. Exchange of topological quanta 

An example of exchange of topological quanta between neighbouring bands is shown in 
Fig. [l^. The figure is an enlargement of the band structure around the X point, showing 
an accidental degeneracy between the 3'ed and 4'th band, which occur about ^ = 0.035. 
Also indicated is the Hall conductance of each band in units of e'^/h, found by numerical 
integration. At the degeneracy it is not possible to define the Hall conductance for the 
individual bands. On the Brillouin zone torus there are two X points Xi, X2, and the bands 
have a first order degeneracy in each. The numerical integration shows that two topological 
quanta are transfered from the lower to the upper band, and this is in full agreement with 
the discussion of Sec. [HI D| . Exchange of topological quanta between bands is a common 



phenomena as ^ is varyed, and this particular example has only been chosen as an illustration 
of the general phenomena. 
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3. Transverse conductivity in a disordered vortex phase 

When the vortices come from a thin film of superconducting material, which have many 
impurities and crystal lattice defects acting as pinning centers for the vortices, the distribu- 
tion of vortices will be disordered rather than forming a regular Abrikosov lattice. We expect 
the effect of the disorder to be, to wash out the distinctive features of the band structure, i.e. 
to average out the characteristic fluctuations in a^y, leaving a smooth curve in Fig. |T^ with 
a characteristic dimensionless proportionality constant s{C,), in the form axy = (e^//i)s(^)z/, 
where z/ is the number of electrons in the magnetic unit cell u = UeO? = ne(f)o/B (the filling 
fraction). With this conjecture we can estimate the normalized Hall conductivity s(^) by 
making a linear fit to the calculated o'xyii') distribution. In the experimental situation A^ is 
constant, and this makes s(^) a function of the applied magnetic field through the relation- 



ship ^ = As/a = Xs\/B/(j)Q. In Fig. |T^ we have plotted s{B) for a vortex exponential length 
Xg = 80nm. The s(i?)-curve shows essentially that the Hall effect of a dilute distribution of 
vortices is strongly suppressed compared to the Hall effect of a homogeneous magnetic field 
with the same average strength. This is in good qualitative agreement with what is seen in 
the experiments of Geim et al.a. When doing experiments, one is not directly measuring the 
conductivities, but rather the resistivities p^x, Pxy The experiments of Geim et al. cover 
the parameter range from Xp <^ Xs at high 2DEG densities, down to the value Xp/Xg = 0.7 
for the 2DEG with the lowest density experimentally obtainable, where the new phenomena 
begin to occur. Our numerical calculations belongs to the other side of this cross-over where 
Xp ^ A^. The physical picture of this cross-over can be stated as follows. On the high 
density side the magnetic field varies slowly over the size of an electron wave packet for 
electrons at the Fermi energy, with the result that the wave packet more or less behaves as 
a classical particle. On the other side of the crossover Xp ^ A^ the magnetic field varies 
rapidly over the length-scale of a wave packet, for an electron at the Fermi energy, and this 
introduces new phenomena of an essential quantum character. 
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F. Superlattice potential 

The general picture we have outhned so far of energy bands having dispersion, with the 
dispersion giving rise to a non trivial behavior of the Hall conductivity, is not limited to 
the inhomogeneous magnetic field. The dispersion could have another origin for instance 
a superlattice potentiaS. To illustrate this a series of calculations have been made on a 
2DEG in a homogeneous magnetic field, and a scalar potential which we have taken as a 
square lattice cosine potential. This system have commensurability problems because of the 
two "interfering" length scales, given respectively by the magnetic length Ib = ^/fijeB^ and 
the period of the superlattice potential a. To make things as simple as possible we have 
fixed the period of the cosine potential a, and the magnetic field strength 5, and only varied 
the amplitude of the cosine potential. Furthermore the flux density of the magnetic field 
is tuned so that the flux through one unit cell of the cosine potential is exactly one flux 
quantum. The cosine potential is 

U{x,y) = Vo(cos27rx/a + cos27r|//a), (94) 

and the magnetic field B = (^q/o? . The dimensionless parameter controlling the shape of 
the energy band structure is in this case 



hjJr 



(95) 



Examples of the cr^-spectra are shown in Fig. [T^. It is observed that although the spectra 
look different from the vortex lattice spectra, they have the same spiky nature. The spikes 
have the same interpretation as in the vortex lattice system. Local spikes are due to local 
gaps in the spectra. That is when to bands are close to each other for some k vector in the 
Brillouin zone, the result is a pile up of topological charge across the gap, and this gives a 
spike in the (Ji^^-spectra when the Fermi energy is swept across the gap. Global spikes, that 
is, spikes which go all the way up to the diagonal line indicating the Landau limit, are due 
to global gaps in the energy spectrum, combined with the topological quantization. 
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IV. CONCLUSION 

In Sec. H of this paper we have considered the longitudinal and the transverse resistivities 
of a 2DEG in a disordered distribution of flux vortices, within the theoretical framework 
where each scattering event is treated independently, and the electrons are non-interacting. 
The general features observed in experiments are in agreement with the results we have 
outlined, but we do not have perfect quantitative agreement. The radius of the vortices, is 
estimated by Geim to be tq — lOOnm, while we find the best fit between the calculated and 
the measured Hall resistivity for ro — 50nm. 



In Sec. |TT| we have considered a new kind of experiment where a 2DEG is placed in a 
periodic magnetic field varying on a length scale A^, comparable to (or less than) the Fermi 
wavelength Xp of the electrons. In this limit, where it is necessary to include explicitly the 
vector potential in a quantum treatment of the electron motion, we expect the 2DEG to 
exhibit new phenomena. We have presented numerical results for a non-interacting 2DEG 
without impurities showing characteristic spikes of the Hall conductivity versus filling frac- 
tion, which can be understood in terms of local and global energy gaps in the spectrum. 
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APPENDIX A: SOLUTION OF THE BOLTZMANN EQUATION 

The Boltzmann equation, linearized in the external electric field, reads 

~ ^'^ ' ^~de " ^" / 12^ {fp+i'^P+q^P - fpWp^P+q} - -^ • (Al) 
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Here n^ is the density of vortices, and Wk^k' is the scattering probabihty for scattering on a 
single vortex. The transition probabihties Wp^p^g are potentially asymmetrical quantities, 
due to the time reversal breaking magnetic field in the vortices. As argued by B. I. SturmanEj, 
the correct form of the collision integral, even in the absence of detailed balance, is the 
one given in Eq. [Al|. The electron- vortex scattering will be elastic. In order to solve the 



Boltzmann equation we Fourier transform, and write 

oo 

f{k,e)= J2 e~^"Vn(fc) (A2) 

n=— oo 
oo 

w{k,e)= J2 e-*"V(A;), (A3) 

n=— oo 

where 6 is the angle between k and E. The electron- vortex collision integral is diagonal 
when Fourier transformed, and we get the following equation for the n'th component of the 
distribution function 

- evE-—-{Sn,i + 5n,^i} = -na{wo - w„}fn —. (A4) 

(76 Z Timp 

The current is given by 

J (27r)2 TThW i^i t 1 

I Wi — Wo ) 

from which the conductivities can be read off. The resistivities are found by inverting the 
conductivity tensor 

u f" rjf) 

P.. = Ph— / —{l + cose)wie). (A7) 

For a vortex with flux a^o we have kp/a = 27 /rg, with 7 = Ic/tq and the cyclotron radius 

Ic = v/uc. 
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APPENDIX B: INTEGRAL REPRESENTATION OF THE AHARONOV-BOHM 
WAVE FUNCTION, AND SCATTERING CROSS SECTION 

In the standard geometry (Fig. ^ the Aharonov-Bohm wave functionta is 

V> = 5^(-2)l— -I J|„_„|(A;r)e*'"^ (Bl) 

m 

By introducing the integral representation of the Bessel functions, the sum can be effectuated 
with the resuh 



(7t 

7/, = / _p-ikrcosT 



1 _ Qi(T-e) I _ g-i(T+e) 



(B2) 



where the contour 7 can be deformed into the contour shown in Fig. |T8|. The integral can 
now be separated naturally into two terms. The integral along the real axis consists of a 
principal value integral which vanish because the integrand is odd, and a contribution from 
the two poles 

^-ikrcose+iae ^^i^ (B3) 

which is nothing but the incoming part of the wave function. We therefore immediately 
have the following expression for the scattered part of the wave function 

IS o • f^ dt .. -,e*^coshat + cosh(l-a;)t ,^^, 

Jo 27r cos 9 + cosh t 

At this level the splitting ■?/' = ^* + ?/'* is exact. Although it may not be apparent, the cut 
in ip^ for 9 = ±7r is accompanied by a counter cut it ip'^, as is necessary in order for ip to be 
smooth. In the "wave zone" kr ^ 1 it is possible to simplify ip'^ by expanding the integral 
asymptotically 



r = - sin ^^e-^rcose^^o/2GiV^cos9/2)^ ^^^^ 

where G, which is related to the error function, is given by 

/■OO 

G{x) = I dse''\ (B6) 

39 



X 



This expression can be further simphfied when \/2kr cos 0/2 ^ 1, i.e., asymptotically in all 
but the forward direction 

ikr+id/2+i-K/4: ikr+in/A 

jj' = -sm7ra , = -F{9) ^, (B7) 

V2^ffc^cos0/2 v^ ' 

leading to the scattering cross section 

1 ^/^m9 1 sin^ vra ,^ , 

Although this differential cross section is singular in the forward direction, it is completely 
symmetric and can therefore not give rise to a Hall effect if substituted into the Boltzmann 



transport equation, Eqns. |A6| - |A7 . 



APPENDIX C: INTERFERENCE BETWEEN THE INCOMING AND 
SCATTERED PARTS OF THE WAVE FUNCTION 

In this appendix we want to argue that it is not the ordinary scattering cross section 
which should be used in the Boltzmann equation, but rather a redefined "cross section" 
which take into account the interference between the scattered and incoming parts of the 
electron wave function. 

The particle density current, of the total wave function ■?/' 

j = Re[ij*vij] (CI) 

can naturally be separated into the terms 

j = Re [tP^vtPi] + Re[V>>V> +V>>V>4 + Re [ipy^ps] • (C2) 

ji jis js 

In terms of the density current the ordinary scattering cross section is 

^ = hm !^^^M, (C3) 

where the index r refers to the radial component, and jo is the incoming current. As is well 
known the scattering cross section satisfies the optical theorem 
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a = 2J^Im[F(7r)], (C4) 



where F(9) is defined as in Eq. |B7[ The optical theorem follows directly from particle 
number conservation V • j + dp/dt = 0. Inside a conductor there are no collimators, or any 
other devises to separate tp^ from ■?/'* at large distances from the scatterer, and consequently 
the physically correct form of the scattering cross section in this case is 

f = hm '^'^-'^ = hm <^r+Jl')_ (C5) 

The equivalent of the optical theorem for this modified cross section is 

a = 0. (C6) 

The interference contribution rff/j^ is only different from zero in the forward direction 9 = 
±7r, because this is the only direction with constructive interference. For all Hamiltonians 
with the mirror symmetry 9 -^ —9, the difference between a and a will be proportional to 
6{9 — tt), and will never show up in the transport coefficients pxx and p^y, because of the 
geometrical factors: 

/7~ 
d9{l + cos9)^ (C7) 

/rirr 
d9sm9— = 0, (C8) 

where l + cos0 ~ (tt — 0)^/2 for 9 r^ n, thus regularizing the integral. But if the Hamiltonian 
is asymmetric in 9, the interference term can have a contribution which is proportional to 
the derivative of a delta function 6' {9 — n), as we will now show. We do this by directly 
calculating the quantity 

lim -^- / d9 sin 9 f;{r, 9) = A^^(a). (C9) 

r^oo 2najo J_^ 

After a some algebra, and dropping terms which integrate to zero, we have 

/tt* sin TTfy I r 

A^^(a) = / d9 sin 9 cos 9 sm(l/2-a)9 S(V2k^ cos9/2) - C(V2k^ cos9/2) , 

Tia J_^ L J 

(CIO) 
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Where 5", C are the Fresnel sine and cosine integral^. For fcr -^ oo we can evaluate the 
integral asymptotically 

.AT,. . 4 sin 7ra COS ttq; f°° , ^^, , ^, ,^ sin27ra /^,,n 

A^^(a) = / dxx[C{x) - S{x)] = ^- . (Cll) 

vra Jo 27ra 

This is exactly the residual value of (^y in the Aharonov-Bohm limit — due to the lordanskii 
force. We conclude that, in this particular example, using the Boltzmann transport equation 
with the ordinary (quantum mechanical) scattering cross section, leads to the erroneous 
conclusion of a vanishing Hall resistance. In order to obtain the correct result it is necessary 
to redefine the cross section to take into account the interference between the scattered and 
incoming parts of the electronic wave function. 
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FIGURES 
FIG. 1. The scattering geometry for classical scattering on an idealized cylindrical vortex with 

constant magnetic field inside and no field outside. Here the radius of the vortex has been taken 

as the unit of length. 

FIG. 2. Differential cross section (divided by ro), and classical trajectories for four different 
values of the parameter 7 = Ic/r^ = 0.025, 0.40, 1.00, 2.50 

FIG. 3. The Hall efficiency factor Cxy{l)^ ^-^d the resistance efficiency factor Cxx{l) calculated 
from the classical cross section for scattering on a magnetic flux tube. The parameter 7 is given 
by the cyclotron radius divided by the radius of the flux tube 7 = Iq/tq. 

FIG. 4. The geometry of the scattering situation. 

FIG. 5. The efficiency of a dilute distribution of vortices in producing Hall effect, com- 
pared to a homogeneous magnetic field with the same average flux density. Curves show C,xy for 
a = 1/4,1/2,3/4,1. 

FIG. 6. Resistance efficiency of single vortex. Curves show 2TTaC,xx^ for a = 1/4, 1/2, 3/4, 1. 

FIG. 7. These plots of Qxy and 2TTaC,xx for a vortex with a = 10, shows a striking structure of 
resonances at the values of kvQ corresponding to the Landau quantization energies. 

FIG. 8. (A) The unit cell with basis. The large circles indicates the position of the Abrikosov 
vortices, and the small circle indicate the position of the Dirac vortex with the counter flux. The 
micro lattice shown here is 6 x 6, whereas all the numerical results we have presented are obtained 
with a micro lattice of 10 x 10 sites. (B) Four concatenated unit cells, showing the square lattice 
of Abrikosov vortices. 

FIG. 9. The Brillouin zone with symmetry labels. 
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FIG. 10. Schematic energy band crossing, controlled by an outer parameter 7 oc ^ — ^o- 

FIG. 11. Band structures for 2D electrons in a square lattice of Abrikosov vortices with oi = 1/2. 
The different plots show band structures corresponding to various values of the parameter ^ equal 
to the ratio between the exponential length of the magnetic field from a single vortex, and the 
lattice parameter. 

FIG. 12. Calculated Hall conductivity versus filling fraction, for various values of the ratio 
^ = As/a. These calculations are made on the same system as the band structures shown in 
Fig. 11, that is, a square lattice of Abrikosov vortices with a = 1/2. Each of the spectra are 
made by calculating the total Hall conductivity ani^p) and the integrated density of states i^iep), 
for 2000 equidistant values of the Fermi energy ep- The x-axis indicates the integrated density 
of states in units of filled bands. The y-axis indicates the total Hall conductivity in units of the 
conductivity quantum e^//i. The diagonal line in the plots indicate the Hall conductivity in a 
homogeneous magnetic field. 

FIG. 13. The density of Hall effect, or "topological charge", plotted as function of the filling 
fraction, for a = 1/2. 

FIG. 14. Band structures for 2D electrons in a square lattice of vortices carrying one flux 
quantum each, a = 1. The different plots shows band structures corresponding to various values 
of the parameter ^ equal to the ratio between the exponential length of the magnetic field from a 
single vortex, and the lattice parameter. The band structures have been calculated using the basis 
shown in Fig. 8, with the only difference that here the fiux through each of the vortices are (po, 
and the counter Dirac flux is — 20o- 

FIG. 15. Exchange of topological quanta. The figure shows an enlargement of the a = 1/2 band 
structure around the X point, where a degeneracy between the 3'rd and 4'th band occur. (See 
Fig. 11). The parameter 7 appearing in the figure is defined as 7 = .^ — ^O) with ^0 = 0.035. The 
numbers give the Hall conductance of the bands in units of e^//i, found by numerical integration. 
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FIG. 16. The normalized Hall conductivity s{B) for vortices of exponential length As =80nni, 
and Q = 1/2. 

FIG. 17. Calculated Hall conductivity versus filling fraction, for a 2DEG in a homogeneous 
magnetic field, and a square lattice cosine potential, in the special case where the magnetic flux 
density is exactly equal to one flux quantum per unit cell area. The dimensionless parameter 
V indicated in the plots, is equal to the amplitude of the cosine potential divided by the Landau 
energy hujc. The x-axis indicates the integrated density of states in units of filled bands. The y-axis 
indicates the total Hall conductivity in units of the conductivity quantum e^/h. The diagonal lines 
in the plots indicate the Hall conductivity in a homogeneous magnetic field, without any potential. 

FIG. 18. The contour of integration. 
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